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Submitted  to  the  Department  of  Naval  Architecture  and 
Marine  Engineering  on  19  May  1967  in  partial  fulfillment  of  the 
requirements  for  the  Master  of  Science  Degree-  in  Naval  Architec- 
ture and  Marine  Engineering. 

The  object  of  this  thesis  was  to  investigate  the  stability 
of  a  ship  with  automatic  steering  in  a  following  sea  using  the 
nonlinear  equations  of  motion  and  to  compare  these  results  with 
the  results  of  linear  approximations. 

The  solution  of  the  nonlinear  equations  was  obtained  by  a 
technique  of  treating  nonlinear  terms  as  linear  terms  with  non- 
constant  coefficients.  Results  showed  that  when  drift  angles 
become  0.1  radians  or  greater  the  nonlinear  solution  differs 
substantially  from  the  linear  solution.  But  the  nonlinear  solu- 
tion has  a  greater  range  of  stability  in  the  wave  profile  than  the 
linear  solution.  This  leads  to  the  conclusion  than  a  system 
designed  using  linear  theory  provides  sufficient  stability  for 
the  ship  in  a  following  sea. 

The  most  important  effect  which  this  approach  neglected 
was  the  coupling  effect  of  the  roll  or  heel  motion  on  the  yaw 
motion.  Before  a  complete  and  accurate  solution  to  the  steering 
problem  can  be  found  the  coupling  effect  must  be  investigated  and 
some  method  of  expressing  it  analytically  must  be  produced. 
Thesis  Supervisor:     Philip  A.  Mandel 
Title:  Professor  of  Naval  Architecture 
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INTRODUCTION 

Steering  a  ship  in  a  following  sea  and  the  danger  of  the 
ship  broaching  in  a  following  sea  are  problems  which  have  been 
prominent  for  many  years.  Before  the  invention  of  the  autopilot 
the  only  solution  was  a  skilled  helmsman  and  an  experienced 
captain  who  avoided  the  situation  whenever  possible. 

Steering  in  a  following  sea  is  not  difficult  until  the  speed 
of  the  ship  becomes  almost  that  of  the  waves.  Then  the  forces  and 
moments  created  by  the  waves  on  the  ship  are  able  to  act  for  rel- 
atively long  periods.  In  todays  fast  ships  it  is  quite  possible  to 
encounter  this  situation  of  low  frequency  of  encounter  of  ship  to 
wave  crest.  There  are  two  main  solutions  to  the  problem.  The  most 
obvious  and  by  far  the  simplest  to  obtain  is  to  reduce  the  speed 
of  the  ship  which  increases  the  frequency  of  encounter  and  short- 
ens the  period  during  which  the  forces  act.  But  in  todays  world, 
time  is  money  and  to  slow  the  ship  in  following  seas  is  undesir- 
able. 

With  an  autopilot  the  steering  problem  can  be  treated  ana- 
lytically because  the  motion  of  the  rudder  is  a  function  of  the 
yaw  and  yaw  velocity  of  the  ship.  The  first  explanations  of  the 
steering  stability  were  made  by  using  the  linear  terms  of  the 
equations  of  motion.  However,  in  ship  motions  such  as  turning 
maneuvers  the  linear  equations  did  not  give  adequate  analyti- 
cal results.  Research  has  been  done  with  the  nonlinear  equations 
of  motion  and  solutions  have  been  obtained  using  important  non- 
linear terms  of  the  equations. These  solutions  were  found  to  give 
adequate  predictions  of  ship  motion  maneuvers. 

In  a  following  sea  with  small  frequency  of  encounter  the 
motions  may  become  large  and  the  linear  solution  may  be  very 


inaccurate.  Until  now  no  one  has  tried  to  use  the  nonlinear 
equations  to  predict  the  stability  of  a  ship  traveling  in  a 
following  sea. 

The  solution  of  the  nonlinear  equations  of  motion  in  the 
case  of  following  seas  is  possible  only  after  making  several 
simplifying  assumptions.  The  velocity  of  the  ship  is  assumed  to 
be  the  same  as  the  wave  thus  the  frequency  of  encounter  is  zero. 
The  ship  position  in  the  wave  profile  is  constant.  The  forces 
and  moments  created  by  the  waves  on  the  ship  are  determined  by 
the  Froude-Krylov  Hypothesis.  Cross  coupling  terms  of  roll  and 
pitch  motion  are  neglected  in  the  equations.  The  surging  motion 
of  the  ship  is  neglected.  Having  made  these  four  assumptions  it 
is  now  possible  to  solve  the  nonlinear  equations  of  motion  for  a 
ship  in  a  following  sea. 


PROCEDURE 
NOMENCLATURE 

The  motions  of  the  ship  are  based  on  axes  through  the  center 
of  gravity  of  the  ship, Figure  I,  and  the  terms  are  defined  below 
or  in  the  reference  material. 

AC  =  X-'V  Heading  angle  from  direction  of  advance  of  the  waves 
&     -      Angle  of  drift  or  side  slip 

w  =   Coefficient  of  yaw  feedback  in  rudder  control  system 
*  =   Wave  length 

i      =   Yaw  angle  or  heading  error 
r    =   Coefficient  of  yaw  velocity  feedback  in  rudder  control 

system 
*•  =   Horizontal  distance  between  wave  crest  and  the  center  of 

gravity  of  the  ship 
Cfc  =   Block  Coefficient 

y    =   Lateral  force  exerted  by  the  waves  on  the  ship 
n   =   Yawing  moment  exerted  by  waves  on  the  ship 
h   =  Wave  amplitude 
»   =   Rudder  angle 

The  derivatives  and  the  nondimensionalizing  factors  for 
them  are  as  defined  in  references  (6)  ,  (10)  and  (11) . 


Figure  I 
Definitions  of  Position, Distances 
and  Angles 


Instantaneous 
Velocity  V 

Prescribed  Course 


*■  Wave  Direction 


NONLINEAR  EQUATIONS  OF  MOTION 

The  motions  of  the  ship  may  be  represented  by  nonlinear  as 
derived  in  reference  (10).  These  equations  have  been  nondimension- 
alized  in  accordance  with  reference  (11) .  The  coefficients  of  the 
terms  in  the  equations  were  determined  for  the  Mariner  hull  model 
in  reference  (6) .  Using  these  equations  and  including  only  terms 
which  reference  (6)  determined  were  important  or  of  minor  impor- 
tance, the  Y  and  N  equations  are  expressed  as 

Y  Equation 

(m  -  Yv)v  +  (mxG  -  Y^)f  -  Yv  v  -1/6  Yvvvv3  -  (Yy  -  mu)r 

-1/2  YyVVYV2  -  Y<&        -1/6  Ystg    *3  -1/2  Njvv  *v2  =  0 

N  Equation 

(mxG  -  N-)v  +  (Iz  -  Ny-)v  -Nvv  -1/6  Nvvv  v3  -  (N  r  - 

,  ,   C2) 

mxGu)r  -1/2  N^vv  yv2  -  Ug£     -1/6  Hgff  6   3  -1/2  N^vv  Jv2 

=  0 

Substituting  #  =  -v  into  equations  (1)  and  (2)  equivalent 
equations  are  obtained  using  drift  angle  in  lieu  of  lateral 
velocity.  Adding  to  equations  (1)  and  (2)  the  force  and  moment 
excitation  of  the  waves  as  developed  in  Appendix  B 

Y  Equation 

(m  -  Yy)  B     +  (mxG  -  Yf)r  -  Yv  tf  -1/6  Yvvv  #3  -  (Y  f  -  mu)^ 

-1/2  Y^vv  yG2  -  Y6S      -1/6  Ym*3  -1/2  Y $vv   S02  =   FyyS    (3) 

cos27t  sin(2  ffaA  )  -  *flFyy/>*  5  cos2  f  sin(2  ^a/x  ) 

N  Equation 

(mxG  -  N£)  6     +  (Iz  -Nf  )f  -  Nv  G    -1/6  Nvvv  tf3  -(N«p- 

mxGu)r  -1/2  N^vvy/?2  -  N,f  -1/6  Nm  *3  -1/2  N  f  vv  *  &2 

(4) 
=    1/2    Fxx    5    cos2^cos(2^a/>    )    -  V 1/2  JFXX/^ !    cos2'*  cos  ( 

2  Tra/N  ) 

I   =   k/Cb  m'4  tt2    hL/>    2  (5) 
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For  a  ship  with  an  autopilot  having  yaw  and  yaw  velocity 
feedback  control  the  rudder  equation  is 

tf  +  -t  /   =  *  V+  r-*f  (6) 

LINEAR  SOLUTION  OF  EQUATIONS 

The  three  equations  in  three  unknowns  are  sufficient  to  give 

a  solution.  The  difficulty  arises  in  the  fact  that  the  equations 

are  nonlinear  differential  equations  and  an  analytical  solution  is 

not  known.  As  a  first  approximation  to  the  solution  neglect  all 

nonlinear  terms  and  find  the  linear  solution.  This  results  in  the 

following  equations 

Y  Equation 

(m   -   Ya)c     +    (mxG   -   Yf)<p     -    (Y  if  -   mu)f«     -   Yv  B  -   Y*  *      + 

(7) 
JF yyA*S    cos2*?  sin(2  7ra/>k   )   *f  =   Fyy   I    cos2-*   sin(2*ra/*    ) 

N   Equation 
(mxG   -   Ny)<?      +    (I2   -   N^)v     -(Ny-   mxGu)v     -Nv  O    -   Nj  J     + 
1/2^FXXA^5    cos2-*  cos  (2  »r  a  A    )  »f  ■    1/2    Fxx   I    cos2  ^  (8) 

cos  (2  jra/^    ) 

Rudder  Equation 
i+ti     =   jfy  +  ry  (9) 

Using  letter  designations  for  the  coefficients  and  the  values 
for  the  coefficients  as  determined  in  reference  (6)  the  equations 
become 

A<3+Cy  +  Dtf  +  E  y  +  G  $  +  II  Y  =  I 

JC+Ky    +    LO+My    +    Pi+Ry^T  (10) 

6  +   t  6      =    y  y  +  <r  y 

where  the  values  of  letter  coefficients  are  given  in  Table  I. 


Table  I 
Values  of  Coefficients  for  Mariner  Hull 

A  =  (m  -  Y^)     =   .01546  J  =  (mxG  -  N^)     =  .000227 

C  =  (mxG  -  Yy  )  =  -.000086  K  =  (I2  -  N  *>  )     =  .00083 

D  =  -Yv          =   .011604  L  =  -Nv            =  .00264 

E  =  -(Yy~  mu)   =   .00499  >M  =  -(Nf-  mxGu)   =  .00166 

G  =  -Y$          =  -.002779  P  =  -N<f            m    .00138 

H  =  )F      /^  i  cos2^  sin(2wa/A)  R  =  h^Fxx/i^l    cos2-*  cos  (2*-a/>) 

I  =  Fyy  5  cos2-ysin(25"a/>  )  T  =  %  Fxx  I  cos2^  cos  (2n-a/>.  ) 

The  general  solution  of  the  simlutaneous  equations 

(As  +  D)  3C    +    (Cs2  +  Es  +  H)  Vi  +  G  £•  =  0 

(Js   +   L)  gl    +    (Ks2    +  Ms   +   R)  v.    +   P  it   =.  0  (11) 

-    (  r  s   +  *  )  V^  +(ts   +    1)  *b     =0 

where         Qt  -  c,eS(t  +  Gtes£  +  ee^  +  Gfesvt 

^  =  ^est  + 

6    =  ^eSt  + 

can  be  found  by  using  determiants.  The  rrots  of  the  solution  are 

thus  the  stability  roots  of  the  system.  The  roots  can  be  found  by 

finding  values  of  s  for  which  the  determinant  is  equal  to  zero. 

CONTROLS  FIXED  STABILITY  IN  WAVES 

For  controls  fixed  stability  there  is  no  rudder  term  because 

the  rudder  deflection  is  zero.  The  controls  fixed  stability  is 

determined  from  the  following  equation: 

(AK  -  JC)  s3  +  (DK  +  AM  -  LC  -  JE)  s2 

(12) 
+  (DM  +  AR  -  EL  -  JH)  s  +  (DR  -  LH)   =  0 


Substituting  real  values  of  coefficients  from  Table  I  into 
equation  (12)  the  stability  of  the  system  is  determined  by  use  of 
Routh's  Criterion  for  various  positions  of  the  ship  in  the  wave 
profile  and  for  various  wave  heights  and  wave  lengths.  Results  are 
shown  in  Figures  II  and  III. 
STABILITY  WITH  AN  AUTOPILOT  IN  WAVES 

For  Stability  with  an  autopilot  in  waves  the  solution  is 
determined  from  the  following  equation: 

(AK  -  JC)t  S4  +  ( (KD  +  AM  -  JE  -  LC)t  +  (AK  -  JC) )  S3 

+  ((DM  +  RA  -  JH  -  LE)t  +  (KD  +  AM  -  JE  -  LC)  +  (AP  - 

(13) 

JG)  <r  )     s2  +  ((RD  -  LH)t  +  (DM  +  RA  -JH  -  LE)  +  (PD  - 

GL)^  +  (AP  -  JG)  y   )  s  +  (RD  -  LH  +  (PD  -  GL)  X   )  =0 
Substituting  real  values  of  coefficients  from  Table  I  into 
equation  (13)  and  assigning  values  of  X+*~  feedback  coefficients 
the  stability  of  the  system  can  be  found  using  Routh's  Criterion. 
Results  are  shown  in  Figures  IV,  V  and  VI  for  the  stability  of 
the  system  in  the  wave  profile  for  varying  values  of  yaw  feedback 
coefficient, with  various  wave  heights  and  wave  lengths. 

In  computing  the  data  for  Figures  IV,  V  and  VI  it  became 
obvious  that  yaw  velocity  feedback  had  no  effect  on  the  range  of 
stability  in  the  wave  profile.  Some  method  of  showing  the  effect 
of  yaw  velocity  feedback  on  the  magnitude  of  stability  was  need- 
ed. This  was  done  by  assuming  negligible  time  lag  in  the  system. 
This  makes  one  root  very  large  and  negative  thus  effectively  re- 
ducing the  equation  to  a  cubic  for  which  the  solution  of  the 
important  stability  roots  was  obtained.  In  solving  the  cubic  only 
the  value  of  the  largest  stability  root  is  determined.  Reducing 
equation  (13)  to  a  cubic  and  solving, the  effect  of  yaw  velocity 
feedback  on  stability  is  as  shown  in  Figure  VII. 
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NONLINEAR  TERMS 

Thus  far  the  affect  of  yaw  feedback  coefficient  and  yaw 
velocity  feedback  coefficient  has  been  shown  using  the  linear 
equations  of  motion.  Nov;  some  method  of  analytically  expressing 
nonlinear  terms  must  be  developed  in  order  to  shov;  their  effect 
on  the  stability  solution  of  the  problem. 

Going  back  to  equations  (3)  and  (4) ,  the  nonlinear  terms  of 
any  importance  are  found  to  be  the  ty  &  2 ,   a  3  an$    6&2    terms. 
Treating  the  terms  as  linear  with  nonconstant  coefficients  they 
can  be  expressed  as  follows: 

Yvvv  /j3  becomes  (  B  2   Yyvv)  B 
Y£vv  S02     becomes   (  e2   Y^vv)  £ 
etc. 

All  the  nonlinear  terms  can  be  treated  in  the  above  form. 
Putting  these  terms  into  equation  (13)  after  assigning  letter 
coefficients  as  shown  in  Table  II  the  following  determinental 
equation  is  produced: 

(AK   -    JC)t    s4    +    ((KD   +    KD1  O 2    +   AM   +   AQ  B2    -   JE 
-JF82    -    LC    -    L'C<32)t   +    (AK   -    JC))     s3    +    ((DM   +    D'Mfi2 
+    DQC2    +    D'Q  <34    +    RA-    JH    -    LE    -    L'E  £2    -    LF  £  2 
-L*F  G4)£    +     (KD   +    KD'    &2    +    AM   +    AQ  O2    -    JE    -    JF  B 
-LC-L,C02)    +     (AP    +    AP,<?2-JG-JG,^2)&')s2       (14) 
+    ((RD   +    RD'  B2    -    LH    -    L'H  *32)t    +    (DM   +    D'M  S2    +    BQ  82 
+    DQ   Q4    +    RA   -    JH    -    LE-L'E  &2    -    LF  /3 2    -    L'FC4)     +     (AP 
+    AP'32    -    JG    -    JG'C2)    *     +     (PD    +    P'DS2    +    PD'G2    + 
P,D'84-GL-GL,62-G,L<32-G,L'04)^r)     s       +  (  ( RD 

+    RD'  3  2    -    LH    -    L'Hfi  2)     +     (PD   +    P'DO2    +    PD'O2    + 
P'D'134    -    GL    -    GL'G2    -    G'Lfl2    -    G'L'G4)^)       =    0 


Table  II 
Coefficients  of  Nonlinear  Terras  for 
Mariner  Hull 
F  =  -1/2  Y^w   =  -.0765    Q  =  -1/2  N  ^w   =   .0274 
D'=  -Yvvv        =   .0808    L'=  -Nvvv        =  -.0164 

G'=  ~Y6vv  "   -0119    p'=  -N6w        ■   -00489 

STABILITY  WITH  AUTOPILOT  IN  WAVES  INCLUDING  NONLINEAR  TERMS 

Substituting  values  from  Tables  I  and  II  into  equation  (14) 
and  by  choosing  an  instantaneous  value  of  drift  angle,  8,    a 
solution  for  the  stability  can  be  determined  for  that  instant  of 
time.  By  choosing  a  value  of  drift  angle  which  is  felt  to  be  the 
largest  expected  value,  the  solution  for  the  roots  is  obtained. 
Plotting  this  solution  along  with  the  linear  solution  a  range  in 
which  the  true  solution  lies  is  found.  These  results  are  plotted 
in  Figures  VIII  and  IX. 

From  all  of  the  previous  equations  it  is  now  possible  to 
investigate  the  steering  stability  of  any  ship  in  waves  and  to 
determine  what  steering  characteristics  are  necessary  to  have 
directional  stability  in  a  following  sea  at  any  given  sea  con- 
dition. 
STABILITY  FOR  MARINER  HULL  IN  REAL  SEA  CONDITIONS 

In  order  to  obtain  information  necessary  to  design  a  steer- 
ing system  that  will  ensure  directional  stability  at  all  times  the 
effects  of  time  lag,  feedback  coefficients  and  rudder  area  on  the 
stability  must  be  known.  There  may  be  limitations  on  the  values  of 
yaw  feedback  coefficient  and  yaw  velocity  feedback  coefficient 
because  operation  in  head  seas  and  calm  water  must  be  part  of  the 
design.  These  operations  may  require  some  range  of  values  of  V  and 
f .  Time  lag  does  not  appear  to  be  a  critical  parameter.  Present 
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system  time  lags  are  in  the  area  of  0.1  and  this  is  sufficiently 
small  for  stability  in  following  seas. 

The  first  consideration  is  that  the  rudder  area  must  be  suf- 
ficient to  counteract  the  moment  created  by  the  waves.  For  the 
Mariner  hull  a  check  for  sufficient  rudder  area  is  made  by  equa- 
ting the  moment  created  by  the  waves  to  the  moment  which  can  be 
applied  by  the  rudder.  Figure  X  shows  the  angle  needed  to  over- 
come moment  of  waves  alone.  The  present  Mariner  rudder  appears  to 
be  large  enough  to  counteract  wave  moment. 

In  all  of  the  developments  using  the  Mariner  hull  the  assump- 
tion is  made  that  the  ship  is  somehow  able  to  maintain  the  same 
speed  as  the  waves.  Thus  the  problem  as  developed  represents  a 
more  dangerous  condition  than  is  actually  present. 

The  next  step  is  to  find  the  value  of  yaw  angle  feedback 
coefficient  necessary  to  make  the  system  stable  in  design  sea 
conditions.  Using  equation  (13)  and  the  values  of  Tables  I  and  II 
the  feedback  coefficient  necessary  for  stability  is  found.  Re- 
sults are  shown  in  Figure  XI. 

Another  question  that  arises  is  how  would  changing  the  rud- 
der area  change  the  stability  feedback  coefficients.  This  is  done 
by  using  equation  (13)  and  varying  the  rudder  force  and  moment 
derivatives.  Using  the  rudder  coefficient  changes  of  reference  (9) 
as  a  basis  a  direct  percentage  change  in  Mariner  coefficients  is 
made  in  order  to  obtain  an  estimate  of  hov;  Mariner  coefficients 
would  change  with  the  same  rudder  area  changes.   The  effect  of 
rudder  area  changes  on  angle  to  overcome  wave  moment  are  shown  in 
Figure  XII.  The  effects  of  rudder  area  changes  on  yaw  angle  feed- 
back coefficient  are  shown  in  Figure  XIII. 

From  the  results  of  these  investigations  in  following  seas 
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and  from  investigations  in  head  seas  and  calm  water  maneuvers  the 
best  combination  of  rudder  area,  time  lag,  and  steering  system 
feedback  coefficients  can  be  chosen. 
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RESULTS 

The  following  figures  show  how  the  directional  stability  of 
a  ship  is  affected  by  the  presence  of  waves  when  running  in  a 
following  sea. 

Figures  II  and  III  show  how  the  range  of  instability  in  the 
wave  profile  for  the  ship  with  no  steering  control  for  varying 
wave  height  and  wave  length. 

Figures  IV,  V  and  VI  show  how  the  range  of  instability  can  be 
decreased  and  even  eliminated  by  the  addition  of  automatic  steer- 
ing controls  with  yaw  angle  feedback  control. 

Figure  VII  shows  how  the  real  and  imaginary  part  of  the 
largest  stability  root  varies  in  the  wave  profile  for  a  system 
which  is  stable  all  the  time  for  two  values  of  yaw  velocity  feed- 
back coefficient. 

Figures  VIII  and  IX  show  how  the  nonlinear  terms  in  the 
equations  of  motion  affect  the  largest  stability  root  and  the 
range  of  stability  in  the  wave  profile. 

Figure  X  shows  the  rudder  angle  necessary  to  counteract  the 
constant  wave  moment  acting  on  the  Mariner  hull. 

Figure  XI  shows  the  minimum  required  yaw  angle  feedback 
coefficient  ncessary  for  a  stable  system  as  a  function  of  the 
sea  condition. 

Figure  XII  and  XIII  show  how  the  rudder  angle  reuired  to 
counteract  wave  moment  and  the  required  yaw  angle  feedback  coef- 
ficient are  changed  when  rudder  area  is  changed. 

Figure  XIV  shows  how  the  minimum  required  yaw  angle  feedback 
coefficient  varies  as  the  size  of  the  ship  changes  for  operation 
in  the  same  sea  conditions. 
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Figure  II 
Stability  of  Controls  Fixed  Ship  in  Waves 
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Figure   III 
Stability   of   Controls   Fixed   Ship    in  Waves 
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Figure   IV 
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Stability  of  Autopiloted  Ship  with  Yaw  Angle 
Feedback  Control 
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Figure  V 

Stability  of  Autopiloted  Ship  with  Yaw  Angle 
Feedback  Control 
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Figure  VI 
Stability  of  Autopiloted  Ship  with  Yaw  Angle 
Feedback  Control 
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Figure  VII 

Real  and  Imaginary  Parts  of  Largest  Stability 

Root 
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Figure  VIII 
Real  Part  of  the  Largest  Stability  Root 
r=  1  y=  1  >/L  =  2   h/A  =  1/40  ^  =  .1 
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Figure  IX 
Effect  of  Nonlinear  Terms  on  Range  of  Stability 
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Figure  X 
Rudder  Angle  Needed  to  Overcome  Wave  Moment 
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Figure  XI 
Yaw  Angle  Feedback  Coefficient  Necessary  for 
for  Stability  Using  Linear  Theory 
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Figure  XII 
Rudder  Angle  Required  to  Counteract  Wave  Moment 
for  Different  Rudder  Areas 
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Figure  XIII 
Required  Yaw  Angle  Feedback  Coefficient  for 
Different 
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Figure  XIV 
Yaw  Angle  Feedback  Coefficient  Required  for 
Stability  for  Various  Ship  Lengths  Using 
Linear  Theory 
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DISCUSSION  OF  RESULTS 
NONLINEAR  EQUATIONS 

The  solution  of  the  nonlinear  equations  presented  an 
impossible  task.  To  date  the  only  solutions  to  the  equations  have 
been  by  numerical  analysis  which  was  designed  to  predict  the 
motion  of  the  ship  during  standard  maneuvers  in  calm  water.  To 
obtain  a  numerical  solution  for  steering  in  following  seas  was 
entirely  beyond  my  comprehension. 

The  only  path  left  open  was  to  try  and  find  some  method  of 
obtaining  an  approximate  analytical  solution.  Previous  analytical 
solutions  had  been  obtained  using  only  the  linear  terms  of  the 
equations.  Investigation  of  the  important  nonlinear  terms  resulted 
in  the  discovery  that  they  were  all  a  function  of  drift  angle 
squared  or  drift  angle  cubed.  By  treating  the  drift  angle  squared 
as  a  non  constant  coefficient  equations  evolved  as  linear  with 
non  constant  coefficients.  These  equations  were  solved  by  choos- 
ing a  maximum  expected  drift  angle  and  putting  it  into  the  equa- 
tions as  a  constant.  This  gave  a  stability  solution  for  an 
instantaneous  value  of  drift  angle.  The  maximum  expected  drift 
angle  solution  is  then  the  solution  for   largest  ship  motions. 

Comparing  the  linear  solution  and  nonlinear  solution  of  the 
largest  stability  root, Figure  VIII,  the  nonlinear  solution  shows 
less  stability  than  the  linear  solution.   Figure  IX  shows  that 
the  ship  has  a  greater  range  of  stability  in  the  wave  profile 
than  predicted  by  the  linear  solution.  These  results  correlate 
with  the  nonlinear  solutions  of  the  standard  maneuvers.  In  the 
standard  maneuvers  of  reference  (9)  the  linear  solution  predicts 
less  stability  and  thus  better  turning  characteristics  than  is 
actually  present.  The  nonlinear  terms  represent  forces  and  moments 
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which  are  always  opposing  the  yawing  motion.  As  a  result  the 
range  of  stability  in  the  wave  profile  is  actually  greater  than 
the  linear  prediction  but  the  oscillating  motion  of  the  system 
does  not  decay  as  rapidly  as  linear  theory  predicts. 

The  actual  stability  of  the  system  lies  somewhere  in 
between  the  linear  prediction  and  the  nonlinear  prediction 
because  the  drift  angle  is  always  changing  as  the  ship  moves 
through  the  waves.  The  graphs  are  plotted  for  a  drift  angle  of 
six  degrees.  In  the  turning  circle  maneuver  the  maximum  drift 
angle  v/as  nine  degrees  and  for  the  zig  zag  maneuver  the  drift 
angle  reached  seven  degrees. 

In  designing  a  steering  system  the  use  of  the  linear  theory 
will  provide  more  than  adequate  control  and  so  the  design  is 
sufficient  for  all  conditions. 
DESIGN  OF  STEERING  CONTROLS 

In  the  design  of  an  automatic  steering  system  the  rudder 
size,  steering  equipment  size  and  the  feedback  to  the  control 
circuits  are  the  parameters  which  must  be  determined.  It  is  very 
important  to  be  able  to  see  what  effect  these  parameters  have  on 
the  steering  problem  in  a  following  sea. 

The  rudder  size  and  the  steering  gear  size  are  initially 
determined  by  the  requirements  set  for  the  ship  to  make  the  stan- 
dard maneuvers  in  calm  water. 

With  this  rudder  size  Figures  II  and  III  show  the  stability 
region  for  the  ship  if  it  has  no  controls  and  the  rudder  remains 
fixed  at  zero  deflection.  The  ship  is  unstable  on  the  negative 
slope  of  the  wave  and  stable  on  the  positive  slope.  This  stability 
region  is  unchanged  by  wave  height  but  does  shift  with  varying 
wave  length  because  the  wave  slope  changes. 
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Adding  automatic  controls  to  the  ship  with  its  present 
rudder-,  Figures  IV,  V  and  VI  show  how  the  value  of  yaw  angle 
feedback  coefficient  required  for  stability  varies  for  different 
types  of  following  seas.  The  value  of  yaw  angle  feedback  coeffi- 
cient needed  for  stability  is  very  dependent  upon  the  size  and 
length  of  the  waves.  As  wave  height  increases  the  more  yaw  angle 
feedback  is  needed  to  make  the  system  stable.  The  important  wave 
lengths  are  in  the  tange  1.0  to  8.0  of  ship  length.  Shorter  than 
this  and  the  wave  has  little  affect  on  steering.  Longer  waves 
require  the  same  or  less  feedback  due  to  decreasing  slope  of  the 
wave. 

Yaw  angle  feedback  is  the  controling  feedback  for  stability. 
The  effect  of  adding  yaw  velocity  feedback  to  the  system  is  shown 
in  Figure  VII.  At  zero  frequency  of  encounter  with  the  waves  it 
has  an  unstabilizing  affect.  But  this  is  for  our  own  hypothetical 
situation.  When  the  frequency  of  encounter  is  not  zero  but  has  a 
value  of  0.5  radians/sec.  or  larger  then  L.J.  Rydill , reference  (8) 
has  shown  that  velocity  feedback  helps  stability.  So  in  an  actual 
following  sea  situation  a  value  of  yaw  angle  velocity  feedback  is 
important  for  good  ship  control. 

Taking  the  Mariner  hull  in  a  following  sea  situation  the 
effect  of  wave  forces  on  stability  can  be  shown.  The  effect  of 
various  sea  states  on  a  500  foot  Mariner  is  shown  in  Figures  X 
and  XI.  Figure  XI  shows  the  minimum  value  of  yaw  angle  feedback 
coefficient  for  operation  in  different  sea  states.  Figures  XII 
and  XIII  show  how  the  minimum  value  of  yaw  angle  feedback  coef- 
ficient is  changed  by  changing  the  rudder  area.  Figure  XIV  show 
how  the  required  yaw  angle  feedback  coefficient  differs  for  a 
200  foot  and  a  100  foot  Mariner  as  compared  to  the  500  foot  ship. 
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Having  obtained  the  data  for  required  yaw  angle  feedback 
coefficient  with  rudder  area  variation, this  must  be  added  to  data 
generated  for  steering  in  head  seas  and  in  standard  maneuvers  to 
give  a  true  picture  of  the  limitations  of  the  steering  parameters 
in  the  different  modes  of  operation.  With  this  data  one  should  be 
able  to  design  the  best  possible  steering  system. 

The  results  for  the  real  ship  case  are  highly  dependent  on 
available  wave  spectrum  data.  The  data  used  here  was  for  fully 
developed  ocean  waves  as  defined  by  the  U.S.  Hydrographic  Office. 
In  coastal  waters  it  is  very  possible  that  the  wave  height  to  wave 
length  ratio  will  be  much  larger  for  given  wave  lengths  due  to 
shallow  water  effects  on  the  waves.  This  fact  would  require  much 
larger  values  of  yaw  angle  feedback  coefficient  to  stabilize  the 
ship.  In  designing  a  steering  system   it  is  most  important  that 
accurate  data  on  operating  sea  conditions  be  used. 
LIMITATIONS  OF  DESIGN  PROCEDURE 

There  are  several  conditions  and  effects  which  may  be  of 
importance,  that  are  not  accounted  for  in  this  investigation.  The 
most  important  is  the  cross  coupling  effect  of  the  roll  motion  on 
the  yawing  motion.  No  one  has  yet  been  able  to  express  analytical- 
ly the  effect  of  rolling  motions  on  the  steering  problem.  It  is 
easy  to  see  that  when  the  ship  rolls  or  heels  during  yawing  motion 
that  it  becomes  unsymetrical  about  the  reference  centerline  and 
side  forces  and  yawing  moments  are  created.  The  larger  the  roll 
or  heel  angle  the  greater  the  forces  and  moments.  This  effect  is 
greatly  dependent  on  the  static  stability  of  the  transverse  motion, 
Investigation  by  the  Germans,  reference  (12) ,in  capsizing  model 
tests  shows  that  as  the  static  stability  of  the  ship  is  reduced 
it  becomes  less  stable  in  yawing  motions  until  it  reaches  a  point 
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where  the  model  broaches  and  eventually  capsizes.  From  their 
investigation  the  importance  of  the  yaw-roll  coupling  is  shown. 
Until  someone  develops  a  method  to  express  this  effect  a  complete 
solution  to  the  steering  problem  is  not  really  known. 

A  condition  which  has  not  been  investigated  is  the  condition 
where  the  frequency  of  encounter  with  the  waves  is  such  that  the 
wave  moments  are  exerted  when  the  ship  is  yawing  away  from  the 
prescribed  course  such  that  the  motion  appears  as  a  spring  mass 
system  being  excited  at  half  frequencies  of  its  natural  frequency. 
It  seems  possible  that  the  ship  yawing  motion  may  somehow  be 
expressed  as  a  simple  spring-mass-damper  system  with  some  natural 
frequency.  Then  the  system  could  be  checked  at  the  half  frequen- 
cies to  see  if  the  motion  is  stable  or  not.  By  changing  the 
ammount  of  feedback  the  natural  frequency  would  be  changed  and 
data  could  be  obtained  which  may  place  further  limitations  on 
the  values  of  the  feedback  coefficients. 

Another  effect  which  has  been  neglected  is  that  of  the 
wind  force  and  moment  on  the  ship.  There  may,  in  a  sufficiently 
strong  wind,  be  forces  large  enough  to  affect  the  steering 
problem. 

The  effect  of  the  surging  motion  of  the  ship  in  a  following 
sea  is  another  condition  which  may  have  important  adverse  affects 
on  the  steering  problem.  For  this  investigation  the  surging 
motion  was  neglected  in  order  to  simplify  the  mathrnatics  of  the 
problem.  The  surging  may  cause  sufficient  changes  in  forces  and 
moments  so  as  to  cause  instability. 


26 


CONCLUSIONS 

1.  The  nonlinear  equations  of  motion  will  provide  a  more 
accurate  solution  v/hen  operating  in  a  following  sea  because  drift 
angles  can  be  expected  to  be  large. 

2.  The  nonlinear  solutions  show  that  when  designing  a  steering 
system  use  of  the  linear  analysis  will  provide  more  stability  than 
is  actually  required  for  operation  in  a  following  sea. 

3.  At  very  low  frequencies  of  encounter  the  yaw  angle  feedback 
in  the  only  control  needed  for  stability.  Yaw  velocity  feedback 
at  low  frequencies  has  an  unstabilizing  effect  but  is  important 
for  good  stability  when  the  frequency  of  encounter  is  above  ■ 
0.5  radians/  sec. 

4.  Values  of  feedback  coefficients  can  be  determined  in  order 
to  make  the  system  stable  when  operating  in  following  seas. 

5.  Values  of  feedback  coefficient  are  highly  dependent  on  the 
wave  spectrum  in  which  the  ship  is  designed  to  operate. 

6.  The  roll  motion  influence  on  the  yaw  motion  can  be  very 
large  for  ships  with  low  tranverse  stability. 

7.  The  point  where  wave  excitation  is  at  half  frequency  of  yaw 
motion  natural  frequency  may  be  as  dangerous  a  position  as  the 
zero  frequency  of  encounter  point. 
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RE  COMMEN DAT I ON  S 
There  are  three  areas  of  the  problem  which  need  further 
investigation. 

1.  Explanation  of  the  Roll-Yaw  coupling  effect  which  will 
produce  the  greatest  changes  in  the  yaw  motion  problem. 

2.  Investigation  of  half  frequency  excitation  which  may  also 
produce  some  startling  results. 

3.  Investigation  of  the  surging  motion  effects  which  may  cause 
slight  changes  in  the  problem. 

It  may  also  be  possible  to  somehow  produce  the  ship  motions 
using  an  analog  computer.  If  this  can  be  done  then  the  nonlinear 
solutions  can  be  produced  exactly  and  the  answer  to  the  half 
frequency  problem  could  be  solved  at  the  same  time. 
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APPENDIX  B 
LATERAL  FORCE  AND  MOMENT  EXCITATION  FROM  WAVES 

The  moments  and  forces  acting  on  the  ship  due  to  the  waves 
are  based  on  the  Froude-Krylov  Hypothesis  which  states  that: 

a.  The  presence  of  the  ship  does  not  disturb  the  pressure 
distribution  in  the  waves. 

b.  The  pressure  distribution  on  the  ship  is  determined  by 
the  hydrostatic  pressure  due  to  the  position  relative 
to  the  wave  crest  with  a  correction  in  order  to  take 
into  account  the  decreasing  orbital  diameters  of  the 
water  particles  with  increasing  water  depth. 

Assuming  a  symetric  ship  with  vertical  frames  at  the  water- 
line  and  following  the  development  in  reference  (1)  the  lateral 
force ,y,  and  the  yawing  moment,  n,  due  to  wave  pressure  are: 

y  -  F  y  •  LBT  •  /»g    2?7h/>  sin(2waA  -  we) 

Having  assumed  the  frequency  of  encounter  is  zero.  we  =  0 

y  =  Fyy  •  LBT  •  /og  2  r?  h~A  sin  (2  77  aA  ) 

n  =  Fxx  *  hi*    '  LBT  •  fiq    2  77-hA  cos  (2^  a/*  ) 

The  coefficients  Fyy  and  Fxx  are:. 

F   =  2/  LB   sin  *    £fw(x)  cos(2"^x/>  )  dx 

--/«. 
Fxx  =  4/L2B   sin  /**v^fw(x)  sin(2^xA  )  dx 

Because   is  not  the  true  angle  between  wave  direction  and 

the  centerline  plane  of  the  ship  the  equations  must  be  modified  as 

follows : 

True  angle  <*  =  *y-  f 

y  -  (Fyy  +  c.Fyy)  LBT  •  pq    2  tt  h A  sin(2^aA  ) 

n  =  (Fxx  +  ^Fxx^  %L   '*  LBT  *  F^   2j»  h/A  cos(2ira/>  ) 
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aFxx  ■   -f  «W** 
y    =    (Fyy   -  f  iFyy/i^  )    LBT    •  />g    2  tt  h/>     sin(2n-a/>    ) 

*    =    *FXX    "V^pxx/J*)     *5L    •    LBT    •    />g    2  7rhA     cos(2  7ra/^     ) 
Non   dimensionalizing    the   equations. 

y»    =  Y/h   u2    S 

n1  =  n/H   u2  SL 
Assuming  that  the  ship  velocity  in  the  direction  of  the 
waves  is  the  same  as  the  wave  velocity 


wave  speed    c  -  V  g  A  /2n 
ship  speed    u  =  c/  cos  // 

u2  =  gA  /2tt  cos2^ 
y'  =  (Fyy  -  r   JFyy/^  )  I  cos2/*  sin(27r  a/A  ) 
**'  "  (Fxx  "  *  iFxx/i//  )  J5  i   cos2,*  cos(27ra/A  ) 

5  =  k/Cb  m1  4  t?  2  hL/p.  2 
k  is  the  factor  used  to  take  into  account  the  influence  of  the 
decreasing  orbital  velocities  as  shown  in  reference  (1) . 

On  the  following  page  are  the  graphs  showing  the  values  of 
Fyy  and  Fxx  coefficients  for  the  Mariner  hull. 
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Figure  XV 
Coefficient  for  Sway  Exciting  Force 


Figure  XVI 
Coefficient  for  Yaw  Exciting  Moment 


W.   (OsjfC'i) 


34 


I 


